Prepared in the fall of 1994 for the proceedings of the 
19th International Nathiagali Summer College on Physics and Contemporary Needs, Nathiagali 1994 

edited by S. A. Ahmad and S. M. Farooqi 
(should have been published by Pak Book Company, but did never appear in print) 



ELEMENTS OF MICROMASER PHYSICS 



Berthold-Georg Englert 

Sektion Physik, Universitdt Miinchen 
Am Coulomhwall 1, D-85748 Garching^ 



Abstract 

The elements of micromaser physics are reviewed in a tutorial way. 
The emphasis is on the basic theoretical concepts, not on technical 
details or experimental subtleties. After a brief treatment of the 
atom-photon interaction according to the Jaynes-Cummings model, 
the master equation that governs the dynamics of the one-atom maser 
is derived. Then the more important properties of the steady states 
of the one-atom maser are discussed, including the trapped states. 
The approximations, upon which the standard theoretical model of 
the one-atom maser is based, are exhibited. The methods by which 
one calculates statistical properties of the emerging atoms are hinted 
at. 



This is a tutorial review of micromaser physics, mainly from the point of 
view of a theoretician. The audience aimed at consists of non-experts who 
want to get an idea about this field of research without perhaps having the 
time to go into the finer details. More technical reviews are, of course, 
available in the published literature. I would like to draw special attention 
to the recent articles of the Garching group and of the Paris group [Q], 
which are both contained in an up-to-date collection of essays on cavity 
quantum electrodynamics In these articles, the reader will find very many 
references to papers on various subtle aspects of micromaser physics that are 
well beyond the scope of this tutorial review. There are also textbooks on 
quantum optics that devote some space to a treatment of the micromaser. 
Some more recent work that is not included in these reviews or textbooks 
will be cited in the text. 



"Also at Max-Planck-Institut fur Quantenoptik, Hans-Kopfermann-Strafie 1, D-85748 
Garching. 



1 




velocity 
selector 



laser 
excitation 



resonator with 
TE mode 



state selective 
ionization 



upper 
level 



lower 
level 



FIG. 1: Schematic setup of micromaser experiments. 

Experimental setup and other introductory remarks 

The most important parts of the setup of a micromaser experiment are 
sketched in Fig. 0. A beam of neutral alkaline atoms (rubidium is the most 
popular species) is run through a resonator, which has a very large quality 
factor to ensure a long photon lifetime. The high quality is achieved by su- 
perconducting surfaces, for which purpose the resonator is made of niobium 
and cooled below the critical temperature of 4.2 K. This low temperature 
has the additional benefit that the number of thermal photons is very small 
(remember the Planck radiation) and thus the thermal noise is reduced. 

The resonator is tuned such that it possesses a mode in the microwave 
regime (of TE type in Fig. |1|) which is resonant with an atomic transition 
between two Rydberg levels. It is this privileged mode that we shall be 
interested in exclusively. The frequencies of all other resonator modes are 
very far away on the scale set by the linewidth of the privileged mode, which 
is extremely narrow owing to the high resonator quality. 

We refer to said atomic transition between the two distinguished Ryd- 
berg states as the masing transition. Prior to entering the resonator the 
atoms are excited to (one of) these Rydberg states. Transitions between 
Rydberg states with principal quantum numbers in the 50s to 70s are in 
the convenient microwave range, where resonators of ultrahigh quality are 
available. In addition, Rydberg atoms are very large (diameters of a few 
thousand Angstrom units are typical) and therefore possess enormous dipole 
moments, so that a strong coupling of the atom to the photons in the priv- 
ileged mode is achieved. It is then possible that the atoms emit microwave 
photons with a large probability. Since the photons are stored for a rather 
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long time in the resonator — lifetimes of a considerable fraction of a second 
have been realized experimentally — even a faint beam of atoms can pump 
the resonator effectively. In this way, a maser is operated in which single 
atoms traversing the resonator provide for an efficient pump. One is then 
dealing with a microscopic maser indeed, a micromaser . 

The excitation of the atoms to the Rydberg states of interest is typically 
achieved by a laser beam which is crossed by the atoms cross just before they 
enter the resonator. It is advantageous to select a narrow velocity group in 
the atom beam because then all atoms interact for the same well-controlled 
time with the resonator photons. A velocity selection can be performed with 
the aid of a set of Fizeau wheels, as indicated in Fig. |l], or by other means 
(see the more technical review articles). 

The high quality of the resonator is quite essential for all micromaser 
experiments. Therefore one cannot perform any direct measurements on 
the stored photon field. The only information available is the state of the 
emerging atoms. Ideally they can solely end up in either one of the states of 
the masing transition. The atoms are probed by state selective ionization in 
inhomogeneous electric fields, symbolically indicated in Fig. |l| by two pairs 
of nonparallel condensor plates. At the first stage only atoms in the upper 
masing level are ionized, in the second stage only those in the lower level. 
The stripped-off electrons are detected and thus the final state of the atom 
is determined. 

The levels of the ^^Rb isotope that are most relevant for many of the 
micromaser experiments in Garching are depicted in Fig. 0. As indicated, the 
masing levels will be called a and (3 throughout. In this figure, the excitation 
from the ground state to the state a is shown as a one-step process. That 
can be done by a suitably tuned, frequency doubled, ring dye laser. But 
different excitation schemes, that involve intermediate levels, have also been 
used successfully. 

The heart of the experiment is the interaction of the excited atom with 
the resonator photons. To a reasonably good approximation only the two 
Rydberg levels of the masing transition are significant for the dynamics dur- 
ing the passage of the atom through the resonator. We shall, therefore, 
simplify matters and think of the masing atoms as two-level systems. The 
state vectors \—) = \a) and \^) = respectively, are then used to denote 
the excited and the deexcited state of such a two-level atom. 
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FIG. 2: Levels of *^Rb relevant for some of the Garching micromaser experiments. 
The laser excitation from the ground state to the upper masing level a is indicated 
by the t arrow, the masing transition by the ]); arrow. The level spacing is not 
drawn to scale. 

Atom-photon interaction (Jaynes-Cummings model) 

A symbolic representation of the passage of a two-level atom through the 
resonator is shown in Fig. ^. Initially the atom is in state a and there are 
supposedly n photons in the resonator, so that the initial state vector of the 
combined atom-field system is \a,n) = |— ,Ti). After the interaction we may 
have either one of two final state vectors — \a,n) if no photon has been 
emitted, or |/3,n + l) = if a photon has been emitted — or a 

superposition of both. 

For a calculation of the actual final state we use the standard Jaynes- 
Cummings model^ for the theoretical description of the atom-field dynamics. 
The fundamental variables of the privileged field mode are the photon ladder 
operators a and with their well known properties, of which 



are perhaps the most important ones. Likewise we employ ladder operators 
a and cr^ for the two-level atoms. Their defining properties are 



[a, a^] = 1 , a^a \n) 
a) \n) 



n) n 



n+l) Vn+ 1 



(1) 





(2) 



^The name derives from a 1963 paper j^]. 
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final state vectors 
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\P,n+l) = \-,n + l) 



FIG. 3: Before and after the passage of a two-level atom that is initially excited 
through a resonator with n photons initially. 

and 

a\(3) = or a\^) = 0, 
a^a) = or a^:^) = , 



with the consequences 
as well as 



= , £7^2 = (4) 

(5) 



l«) ("I = I-) (-1 : 

The latter are useful ways of expressing the projectors to the atomic states. 

In the Jaynes-Cummings model the dynamics is generated by the Hamil- 
ton operator 

H = hujo^a + hVla"^ a — hg{t){a) a + a'^ a) . (6) 
photon atom interaction 

Here, huj is the energy per photon, Ml is the atomic level spacing (the de- 
excited atom is assigned energy zero by convenient convention), and the so- 
called Rabi frequency g{t) measures the strength of the interaction. Inasmuch 
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FIG. 4: Time dependent Rabi frequency g{t) (solid line) and effective Rabi 
frequency g. 

as tlie atom enters and leaves the resonator, this Rabi frequency possesses 
a natural time dependence. It vanishes as long as the atom is outside the 
resonator, and is large while the atom is inside. To the extent to which the 
atomic center-of-mass motion is classical,! the duration of the interaction is 
given by the ratio L/v of the cavity length L and the atomic velocity v. These 
circumstances are sketched in Fig. ^ where we also indicate the effective Rabi 
frequency g, which is given by the time average of g{t), 

-g = ^jdtg{t). (7) 

To be specific let us report typical numbers for the micromaser experiments in 
Garching that make use of the Rydberg transition of Fig. ^ The resonator is 

^This condition is very well obeyed under the ordinary experimental conditions. Under 
extreme circumstances, however, the quantum nature of the center-of-mass motion may 
become important Q. 
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of cylindrical shape with an almost circular cross section of diameter ~24 mm 
and a length L that is also ~24 mm. The experiments use the photons in the 
TE121 mode. Here one has the frequency values 

oj^VL ^ 27r X 21.5 GHz, 

(8) 

(7 ^ 44 kHz = vr X 14 kHz , 
which correspond to the energy values 

hu = hn = 9xlO-^eV, 



(9) 

hg = 3x lO^^eV. 

Please note how small is the energy per photon^ and how tiny is the interac- 
tion energy. 

The transition operator 

7 = aV + a"^a (10) 
effects atomic transitions along with the emission or absorption of a photon: 



7|-,n) = \zi:,n+l) y/n + 1, ^^^^ 
7|^:,n + l) = \:^, n) \Jn + 1 . 

The eigenvalues of 7 are 7' = 0, ±1 ± -\/2, ±-\/3 . . ., and its eigenstates (so- 
called 'dressed states' of one kind in the jargon of quantum optics) are very 
simple linear combination of the atom-field states that we have been using 
so far, viz. 



|7o) = 1^,0), |7^) = ^(|^,n)±|^,n+l)) 
7l7o) = 0, 7l7^) = |7^)(±v/^^) ■ 

They are mutually orthogonal and normalized to unity. 



;i2) 



if 7' = 7" 

(7'l7") = 5(7,7") = <! * ^ ,^ „ (13) 

if 7 7^7 



^Recall that visible photons have energies of 1-2 eV. 
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and complete, 



(14) 



and constitute a convenient basis for expanding arbitrary state vectors. 
The important identity 
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(15) 



is easily demonstrated with the aid of the algebraic properties of the ladder 
operators stated in Eqs. ([l|)-(|D. We use it to rewrite the Hamilton operator 
(P) into a form, 



H = hu-f^ + h{n - u;)o-V - hg{t)'y , 
which is particularly transparent on resonance [uj = Q): 

H{t) = huj'j^ - hg{t)'y. 



(16) 



(17) 



For the sake of simplicity, we shall confine ourselves to this resonant situation 
for the rest of this tutorial review. 

Since the transition operator is then a constant of motion — the vanishing 
commutator [H(t),'y(t)] = implies 7(t) = 7(^0) — the Schrodinger equation 
obeyed by the eigenstates of 7(t) is very simple. For the left eigenstates, it 
reads 

d 



{^',t\H{t)= hui^-hg{t)i (y,t 



and can be solved immediately. The outcome is 

(7',t| = exp[-iu;7'2(t - to)] exp[i7' / ^t' g{t') \ (7',to| • 



to 



(18) 



(19) 



free evolution 



interaction 



We note that the evolution is naturally split into two contributions: the free 
evolution and what results from the interaction. 

We are particularly interested in the overall change asked for in the con- 
text of Fig. 1^. Then the instants to and t are 'before' and 'after' the inter- 
action, so that the integral in (|19|) equals the effective Rabi frequency of (|^) 
times the classical interaction time L/v, 



I dt' g{t') = gL/v = Lp. 



(20) 
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For obvious reasons, the quantity if here introduced is called the accumulated 
Rabi angle. This single number summarizes the net effect of the atom-field 
interaction. 

We have thus found a first central result of micromaser theory: 



Resonant interaction leads to a change 



(Yl 
lY) 



on top of the free evolution. 



(21) 



As an application let us reconsider the situation of Fig. |^. First, we write the 
initial state vector |:?:, n) as a superposition of eigenvectors of the interaction 
operator 7, 



(22) 



= i(|-,n) + |^:,n + l)) + i(|-,n) - 

= ^,{\lt) + \ln)) ■ 

Second, we make use of ( pT]) to find the change resulting from the interaction, 
1 



interaction 



72 



7+) exp(-2(/}Vn + 1) + 7„ )exp(i(y9Vn + 1) . (23) 



Third, we express the outcome in terms of the atom-field states, 

\—.,n) '"^'^'''^'^^'"^ > I — , Ti) cos(v9-\/n + 1) + | m, n + 1) [—i sin((y9 VrT+T)] , (24) 



which is the desired result. In particular, it states the probabilities with 
which the two possible final atom-field states occur. 



. \ interaction 



I— ,n) with probability coi?{ip\/n + 1) 



with probability sin {}p\/n + 1) . 



(25) 



It is essential to realize that this result is dramatically different from 
the one in the usual spontaneous-and-induced emission process. There the 
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physics is well described by the lowest-order perturbation theory, in which 



the argument >~p\/n + 1 of the sine and cosine function in (|25|) is infinitesimal 
and individual transitions are not happening frequently. In addition, usual 
spontaneous emission involves not only just one mode of the radiation field 
but a continuous range of modes with neighboring frequencies. In summary, 
the photon emission in a micromaser experiment has nothing in common 
with the much more familiar spontaneous emission by atoms in free space. 

We note that there are special situations in which the photon emission 
happens with certainty, 

I . \ interaction i i i \ -r / 1 — T /n 

>\iiL,n^\) if v5Vn + l = 7r/2 

(26) 

(e.g. n = 0, Lp = 71/2) , 
or in which surely no emission takes place, 

I . \ interaction i . \ -c / ; — T 

^1-,^) if (/?V«+l = 7r ^^^^ 

(e.g. = 0, ip = Tx). 



Of course, one could add integer multiples to i^^/n + 1 in both cases without 
changing the final states. Clearly, these special situations are unknown in 
the context of the usual spontaneous emission. In the micromaser the atom 
can first emit and then reabsorb the photon — this is the case of (pTf ) — 
thus undergoing a complete Rabi cycle, or perhaps a few complete cycles. 



One-atom operation 

The beam of atoms in Fig. |l] comes originally from an opening in an oven, has 
been collimated, and perhaps velocity selected, before it is directed through 
the resonator. If no special action is taken, as is the standard experimen- 
tal situation, the atoms in the beam are uncorrelated. In other words, the 
atoms arrive at random. The probability p{t)dt that two successive atoms 
are separated in time by t ■ ■ ■ t + dt is thus given by the Poisson formula^ 

p{t)dt = e-^Vdt , (28) 

*The Poisson formula governs also more mudane phenomena, such as the temporal 
spacing between drops of rain. 
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FIG. 5: An atom gives rise to a one-atom event if the temporal separations from 
its predessor, t+, and from its successor, i_, both exceed the classical interaction 
time L/v, equal to the ratio of the cavity length L and the atomic velocity v. 



where r is the rate at which the atoms arrive. The average temporal spacing 
between atoms is clearly equal to 1/r; the average spatial distance is then 
v/r if all atoms move at the same speed v, a condition that we shall assume 
to be obeyed for the time being. 

If the beam is faint enough, most of the atoms will enter the resonator 
after the preceding atom has left it, and will leave before the succeeding atom 
arrives. This is the situation depicted in Fig. |. The probability that the 
middle atom in this figure gives rise to such a one-atom event is equal to 
the product of the probabilities that both spacings, t- and t+, exceed the 
classical interaction time L/v. Now, the Poisson statistics ( pSj) imply that 
the probability for a spacing greater than a certain specified period T is given 
by 

dtp{t) = e-'^. (29) 



(30) 



Consequently, the asked-for probability of a one-atom event is 
exp(— rL/ti) X exp(— rL/f) = exp(— 2rL/t>) . 



t->L/v 



t+>L/v 



A typical number for L/v is 50 /is, so that at the rather large beam rate 
of r = 1000/s this probability is about exp(— 0.1) = 90%. In other words, 
90% of the atoms would contribute to one-atom events. More typical is 
a beam rate of r = 10/s, for which the probability for one-atom events is 
about exp(— 0.001) = 99.9%. Collective events, in which more than one atom 
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interact with the resonator photons at the same time, are then extremely rare 
and can be justifiably disregarded. We are then dealing with the One-Atom 
Maser, the most intriguing version of the micromaser. 

One-atom maser: Dynamics 

For the theoretical treatment of the dynamics of the photon field in the one- 
atom maser (0AM) we employ the standard master equation approach. It 
deals with a coarse grain time evolution of the photon field. Changes in the 
field that occur on the short time scale, on which the atom-field interaction 
happens, are ignored. In essence, one is only interested in the changes that 
occur in the photon state during the passage of a fair number of atoms. 

Since we are focussing on the state of the photon field, rather than on 
the state of the entire system consisting of very many atoms and the photon 
field, we are treating a partial system. More precisely: we are treating an 
open, driven quantum system. Its state cannot be described in terms of a 
state vector | ) or (one of) its wavefunctions. What one needs is a state 
operator p, which is a positive operator of unit trace. These two properties 
of p ensure that all probabilities ( | p | ) are non-negative and that the sum of 
all probabilities of a complete set of mutually exclusive alternatives is unity. 

The matrix elements of p make up the so-called density matrix; for in- 
stance, when the photon number states {n\ and |m) are used, we get the 
number-state density matrix pnm = {n\ p \m). Just like there are very many 
wavefunctions to one state vector | ), there are also very many density ma- 
trices to one state operator p. It is clearly advantageous to concentrate on 
the (somewhat abstract) state operator in the first place, and to choose a 
convenient matrix representation later when it is really needed. 

For the photon field inside the resonator of a 0AM, the state operator p 
— we shall frequently speak of 'the state p' simply — is a time-dependent 
function of the ladder operators and a: p — pt{a\a). For the benefit 
of more transparent equations, we shall suppress the arguments of p where 
misunderstandings are unlikely. The master equation obeyed by p has the 
general structure 



d 


d 




d 




d 




diP 


+ 

free 


dtp 


+ 

gain 


dtp 
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FIG. 6: Gain in 0AM operation. The two-level atom is excited when it enters 
the resonator; the initial atomic state is |:?:)(:?:| = |a)(a| = a'^a. Prior to the 
interaction the photon field is in the state p{a'^,a). Thus the initial state of the 
complete system is P = p{a^ ,a)a^a. The interaction changes this state. We do 
not care about the final atomic state, but want to know the photon state after the 
interaction. 

The 'free' contribution is, of course, generated by the photon part of the 
Hamilton operator (^, 

d 1 

-g^P = [^^photon,p] = -iuj[a''a,p] . (32) 

free 

The 'gain' contribution to (|3lD is the additional change resulting from the 
interacting atoms traversing the resonator, and the 'loss' contribution orig- 
inates in the coupling of the photons to their surroundings. We shall first 
derive the gain term and then turn to the loss term. 

The gain is determined by considering the situation of Fig. |^. After 
employing the number-state matrix of the initial field state p, 

oo 

p(a\a)= J2 \n')pnm{m\, (33) 

n,m=0 
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to write the complete initial state P in the formQ 

P = ^\—,n) pnm{-,m\ . (34) 

Now we can make use of Eq. (^) to find the effect of the interaction, 

Pnm 



p jnteracUon_^ ^ 1 1?:, n) cos(v5 VrT+T) — | n + 1) i sin((/?y?T~TT) 

n,m 



X 



cos{ip\/m + 1) {—, ml + i sm{ipy/m + 1) (m, m + 1| 



(35) 



after ' 



This is the final state operator of the complete system consisting of one atom 
and the photon field (we remember that the 'free' evolution comes on top). 
We do not care about the atom and therefore we trace over the two-level 
degree of freedom (symbolically: tratom) to obtain the final photon state, 

/ t N interaction , r n 1 

p(a' , a) ^tratom {-Pafterj • (36) 

In more detail this reads 



P= T, \n) pnm {m\ 

n,m 

interaction 



J2 \n) cos{(py/n + I) Pnm cos{(py/m + 1) (m| 



n,m 



+ |n + 1) sm{(py/n + l)pnm sm{ipy/m + 1) (m + 1 1 

^(37) 

Upon recalling that a^(aa^) is a normalized ladder operator, 

a^(aa^)-^/2 \n) = \n + 1) , (38) 
one establishes the identities 



\n) cos{(^y/n + 1) = cos[ip{aa^y^'^]\n) 



I W tsinkMffll \ ^^^^ 

\n + 1) sm{ipy/n + 1) = — /^ uih — l^/ ■ 



^The letter P is a capital p, not a capital p. 
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These are then used to arrive at a second central result of micromaser theory: 



Resonant interaction of atoms entering in 
the — state leads to a change 



P 



cos[(p{aa^y/'^] p cos[(p{aa^y/'^] 



+a 



sm[(p{aa 



[aa 



UTJ2 P 



sin[(f{aa 



tU/21 



(aa 



tW2" 



in the photon state on top of the free evolution. 



(40) 



We add the remark that an analogous result with and a interchanged is 
found for atoms that enter the resonator in the ^: state. The derivation is 
left to the reader ClS db simple exercise. 

The gain term in (0) is now available. It is given by 



d 

diP 



beam rate 



change by one atom 



gam 



r< cos [v? (aa"!" ) ^/^] p cos[(y9(aa"l")^/^] 



(41) 



+a^ 



sin[y(aat)V^] sin[^(aat)V^] 1 
P (aan^^ P\ 



aa' 



This expression is valid for the coarse grain time scale provided that the 
arriving atoms are uncorrelated, as we have assumed in Eq. (p8|), for example. 

We turn to the loss term. Here we have to model the coupling to a thermal 
reservoir. In view of what we have at hand, the simplest model is that of a 
thermal beam of very many, very weakly interacting atoms. This beam is 
supposed to be at a certain temperature, so that the rate of — atoms is 
larger than the rate of the atoms. The ratio of these partial rates is, of 
course, given by the Boltzmann factor. 



ra/rp = exp 



Tiuj ^ 



V 



(42) 



where /cb is the Boltzmann constant and T the temperature. Anticipating 
its convenience, we introduce the parametrization in terms of the positive 
number v whose physical significance will become clear later. 
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For the — atoms in the beam, we have a contribution to the loss term 



which is analogous to the right-hand side of Eq. (41), except that we are 
now dealing with a tiny Rabi angle rather than with the macroscopic angle 
ip. The lowest non-vanishing order in is all that we are going to take into 
account. Thus we find 



d 



loss,ci 



' t 1 t 1 t' 

a' pa — —aa'p — -paa^ 



(43) 



for the — contribution. The m contribution is immediately available too, 
because the remark after (|40|) implies that we have to interchange the ladder 
operators a and a\ 

d 



dt 



loss,/? 



apa' 



-a)ap 



-pa) a 



(44) 



Upon introducing the loss rate A = (r^ — To)0^ and the parameter u of 
we combine these partial answers into 



d 

diP 



loss 



a'^ap — 2apa'^ + pa'^a 



aa^p 



2a^pa + paa^ 



(45) 



The three ingredients of the master equation (|3T| ) are now at hand. 

The physical significance of the rate A and the number u is revealed by 
a quick look at the dynamics of an unpumped resonator, for which 



d 

diP 



loss 



— zu;[a"l'a, p] — |y4(z/ + 1) 



-lAu 



a'^ap — 2apa"l' + pa'^a 



aa^p 



2a^pa + paa 



(46) 



is the equation of motion. This master equation is in fact so simple that 
the general solution for an arbitrary initial state can be written down in a 
number of ways. In the present context that is not so essential. We are 
content with these observations: 

(i) The stationary state of (^6|) is the thermal photon state 



Pth 



a'^a 



1 — exp 



Tioj 



exp 



huj 



a'^a 



(47) 
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which identifies v as the number of thermal photons inside the resonator 
at temperature T. 



(ii) For the evolution of the mean photon field (a) = tr{ap} (where tr 
trphoton for short) we have 



-iu - A/2){a) 



(48) 



so that 



(49) 



This shows that the field decays with the rate A/ 2 towards its station- 
ary null value. 



I^iii) Likewise we find for the mean number of photons, 



d 



dt 



{a)a) = tr < a)a—p > = —A{{a)a) — u) 



dt 



(50) 



and therefore 



At 



(51) 



The mean photon number decays with rate A towards its stationary 
value of z/. 

Thus the meaning of A and u is quite clear. 

The relation between the temperature of the resonator and the number 
of thermal photons deserves a little attention. With the numbers of Eqs. (|^) 
and fP) we find the entries in this table: 



T 



1.5K 
0.5 K 
80 mK 



1.96 
0.15 

2.5 X 10-6 



remarks 



early experiments 
typical values 
achieved 



Please observe that in the microwave domain thermal photons are usually 
present in considerable amounts. This situation is quite different from that 
of the optical domain in which thermal photons are irrelevant. 
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[a] pump gain from a^a — 1 

[b] thermal gain from a^a + 1 

[c] thermal gain from a^a — 1 

[d] pump loss to a^a + 1 

[e] thermal loss to a'I'a — 1 

[f] thermal loss to a'I'a + 1 



FIG. 7: Probability fluxes representing the gain and loss terms in the 0AM master 
equation for a diagonal photon state. The labeling agrees with the one in Eq. (|53|). 

Another remark concerns the connection between the decay constant A 
and the quaUty factor Q of the resonator. Per definition the energy stored 
decreases to 1/e-th of its initial value while the field undergoes Q cycles, 

^-At ^ g-i ujt = 27rQ : A = ^ / Q . (52) 

In the Garching micromaser experiments one has u = 2tt x 21.5 GHz, and 
quality factors of Q = 10^ ■ ■ ■ 10^° have been achieved, so that the photon life 
time 1/A can be as large as a few hundred microseconds. 

One-atom maser: Steady state. General matters 

The 0AM master equation (^) with (|32D , (0), and (|45|) does not contain 
any terms that are sensitive to the phase of the photon field.Q Another way 
of putting this states that the diagonals of the number-state matrix of p are 
dynamically uncoupled. Therefore, a state p which at one time is diagonal in 
this sense, stays diagonal for all times. A diagonal state is represented by a 
state operator p that depends solely on the product a^a, the photon number 
operator, but not on the ladder operators a and individually: p = p{a^a). 
The numbers p{n) = {n\ p{a}a) \n) are the probabilities to find n photons in 
the resonator. 

^hi more technical terms this statement says that (^ij) is invariant under the unitary 
transformation p — > pU with U ~ exp(ixa^a) where x is an arbitrary (constant) phase 
angle. 



i^a + 1 



■ a^a 



■ a^a — 1 
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For such a diagonal state, the master equation ( pTD is particularly trans- 
parent. It constitutes a third central result of micromaser theory: 



Master equation for diagonal states 
of the one-atom maser: 



r sin^ 






-1) 
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+A{u + l){a^a + l)ptia^a 
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[B 
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—r sin^ 
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-A(z/ + 






[E] 








[a)a) 


[F] 



(53) 



The six terms on the right-hand side possess obvious meanings. These are 
stated in the probability flux diagram of Fig. |^. 

In the steady state p*-^^^ of (p3D — which is, of course, also the steady 
state of ( piD — the gains and losses of p{a)a) have to add up to zero. A 
rearrangement of the right-hand-side terms in (|53|), 



r sm 



+Ah'a)ap^^^\a^a - 1) 
-A(z/+ l)atap(ss)(^ta) 



r sm 



P 



(SS)|'„t 



(^(ata + 1)1/2 

Az/(ata+ l)p(ss)(ata) 
-A(z/ + l)(ata+l)p(ss)(ata + l) 



[c] 

N 

[F] 



(54) 



combines them into two groups of three. Now, observe that the replacement 
a^a a^a + 1 turns the [a,C,e] group into the [d,F,b] group. By induction 
we therefore conclude that the value of the [a,c,e] group is the same for 
each value of a^a. But for a^a = it is zero. So we arrive at the two-term 
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FIG. 8: In steady state the probability fluxes obey a detailed balance between 
neighboring photon numbers. The terms are labeled the same way as in Eq. (^) 
and in Fig. 0. 



recurrence relation 

A{iy + l)a^ap^^^\a^a) = (r sin^ [if {a^ay/^] + Aua^ a) p^^^\a^ a - 1) , (55) 



which is equivalent to the three-term relation (p^ but considerably simpler. 

The transition from ( p^ to (^) recognizes the detailed balance in steady 
state, illustrated in Fig. |. It comes about because the photon number ladder 
ends at the a^a = rung. 

The recurrence relation (^) can be solved immediately, and we obtain a 
fourth central result of micromaser theory: 



The steady state of the one-atom maser, 
pumped by resonant — atoms, is given by 



JSS) 



(o)n 



n=l L 



+ 



r/A siia^^ipy/n) 



u+1 u+1 



n 



where the value of p^^^^(O) is determined 
by the normalization of p*^^^-* to unit trace. 



(56) 



We note that this steady state is determined uniquely by three numbers: the 
thermal photon number z/, the accumulated Rabi angle if, and the ratio r/A 
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of the pump and the decay rates. This ratio tells us how many pump atoms 
traverse the resonator (on average) during one photon life time. Therefore, 
we call vjA the effective pump rate.[] 

A simple, but necessary, check of consistency is provided by the consid- 
eration of the r = case, the case of an unpumped resonator. The steady 
state of (|56D is then the thermal state pth of (^Tf ), as it must be. 



One-atom maser: Steady state. Examples 

For fixed values of the effective pump rate rjA and the thermal photon 
number z/, the steady state of the 0AM is traditionally studied as a function 



of the so-called 'pump parameter' Q = Lp^r/A. Why this scaling of the Rabi 
angle appears convenient will become plausible later. 

For a certain 0AM steady state p — for simplicity, we drop the superscript 
of p^^^^ until we start considering the time dependence of p again — the 
fundamental statistical quantities of interest are mean the photon number 
(a^a) = trja'l'ap} and its normalized variance, 

— m — ■ ^ ^ 

This quantity (or its square root) is known under various names in the lit- 
erature,0 such as 'Fano factor,' or 'Mandel's Q parameter,' or 'Fano-Mandel 
measure.' We prefer the first one in this list. The normalization in ( |57| ) is to 
the value of a Poissonian distribution (as in a coherent state of the harmonic 
oscillator or a classical state of the radiation field), 

p{n) = ^(a^a)"exp(-(a^a)) , (58) 

Poisson 

for which ((a^a)^) — {0)0}"^ = {a'' a). With this normalization comes the 
jargon of calling a state p 'superpoissonian' if its Fano factor exceeds unity, 
and 'subpoissonian' if it is less than unity. An example for a superpoissonian 
state is the thermal state of Eq. (|47|); its Fano factor equals z/ + 1. 



^In the micromaser literature the symbol A'ex is frequently used as an abbreviation for 
r/A where the subscript 'ex' is a reminder of 'excited.' We do not adopt this custom. 

^These names derive from a 1947 paper by Fano and a 1979 paper by Mandel ||], 
I believe. 
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FIG. 9: Mean photon number and Fano factor in the 0AM steady state for 
an effective rate of r/A — 200 and a thermal photon number oi u — 0.15, as 
a function of the pump parameter 9 in the range from ^ = to ^ = IOtt. 
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FIG. 10: Fano factor in the 0AM steady state for an effective rate of 
r/A = 200 and a thermal photon number of z/ = 0.15, as a function of the 
pump parameter in the range < ^/vr < 10. 



In Fig. 1^ the dependence of the photon number and the Fano factor on 
the pump parameter 6 is shown for an effective pump rate of r/A = 200 and 
u = 0.15 thermal photons. As 6 increases from ^ = 0, we see an initial rapid 
growth of both the photon number and the Fano factor. Both numbers reach 
a maximum and then decrease toward a minimum near 9 = 2ti. This pattern 
repeats itself roughly with a period of 27r. 

The maximal value of {a) a), reached near = 7r/2, is very close to 200, 
the value of the effective pump rate. Thus, for this pump parameter the 
interaction time is just right to ensure that almost all atoms emit a photon 
into the resonator. At the minimum near 9 = 27r, by contrast, a large fraction 
of the atoms do not emit. The interaction time is longer, so that the atoms 
have a good change to undergo a complete Rabi cycle |— ) \^) |— ) and 
leave the resonator in the excited state. 

The bottom part of the plot of the Fano factor in Fig. ^ is enlarged in 
Fig. |lO[ We observe that there are large 6 intervals in which the Fano factor 
is less than one, thus indicating a subpoissonian state of the photon field. It 
is clear that such states can be realized experimentally, since a very precise 
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control of the parameters r/A, z/, and 9 is not necessary to stay inside the 
appropriate parameter range. 



A look at the histograms in Fig. O might be helpful in understanding 
why the Fano factor covers such a large range of values. For reference each 
histogram is compared with a corresponding Poissonian one with the same 
mean number of photons. Whereas we have an almost poissonian distribution 
for 9 = 7i/2 (at the maximum of the (a^a)), the distribution for 9 = 0.324427r 
is much broader than the Poissonian one, and that for ^ = 27r is much 
narrower. In these cases it is easily understood that the respective Fano 
factors are close to, much larger than, or much smaller than unity. A different 
situation is observed for 9 = 2.09 where, according to Fig. |^, the Fano factor 
is particularly large. In Fig. ^ we recognize the reason for that: there are two 
narrow peaks — each subpoissonianly narrow — separated by many photon 
numbers n. As a consequence, the Fano factor, which is a global property of 
the distribution, is enormous. 

The notion of the 'maser threshold' is also most easily understood in the 
context of such histograms. The thermal state (^) is surely below threshold 
by any definition. Likewise the possonian state (|58|) , which is so typical 
for the output of a laser, is certainly above threshold. A striking difference 
between the two is that the maximum of the poissonian state is at a large 
photon number n, whereas in the thermal state the largest probability is 
the one for zero photons. This invites the following criterion based on the 
relative size of p(0), the probability for no photons, and p(l), that for one 
photon: 



p(0) > p(l) 
p(0)=p(l) 
p(0) < p(l) 
In view of the recurrence relation 



-sin2(^) 



1<J 




below 
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above 



below threshold, 
at threshold, 
above threshold. 

, the explicit 0AM version reads 



threshold. 



(59) 



(60) 



In particular, if the effective rate r/A is large, the threshold value of cp is so 
small that sin^{ip) can be replaced by y?^, and then we have simply 9 = 1 
at the maser threshold. This is one reason why the pump parameter 9 is 
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FIG. 11: Histograms for the photon number distribution in the 0AM steady 
state for an effective rate of r/A — 200 and a thermal photon number of 
u = 0.15. The plots are for (a) 9 = 0.324427r (i.e., </? = 4.129°), (b) 9 = 7r/2 
(i. e.,ip = 6.364°), (c) 9 = 2% (i. e., (p = 25.46°), and (d) 9 = 2.09% (i. e., (p = 
26.60°). The solid line shows the actual distribution. The broken line with 
the shaded area underneath shows the corresponding Poisson distribution 
with the same mean photon number. The respective mean photon numbers 
(a^a) are (a) 30.14, (b) 199.44, (c) 36.25, and (d) 52.72. The Fano factors 
equal (a) 12.56, (b) 1.153, (c) 0.214, and (d) 28.68. 
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FIG. 12: Histograms for the photon number distribution in the 0AM steady 
state for an effective pump rate of r/A = 4 and u = 0.1 thermal photons. 
The plot shows the histograms for three different Rabi angle: = 30° (line 

, at threshold), (p = 25° (line , below threshold), and = 35° 

(line — • — above threshold). 



sometimes preferred over the unsealed Rabi angle (p. These matters are 
illustrated in Fig. ^ where the value r/A = 4 is chosen for the effective rate, 
so that (fi = 30° holds at the threshold. We note that the thermal photon 
number u plays no role in (|60|). 

The threshold concept has its obvious limitations when the photon num- 
ber distribution may have more than one peak, as realized in the 6 = 2.09 



case of Fig. 11. More interesting is the question of how many minima and 
maxima are there and at which n values. The answer is obtained after first 
recognizing this consequence of the recurrence relation (^5]): 



if -J sin {^p)\/n) > n then p{n) > p{n — 1) , 

V (61) 

if -^siv? [ip y/n) < n then p{n) < p{n — 1) . 

In this way the n ranges in which p{n) decreases or increases monotonically 
are determined. At the switch-over values we find the maxima and minima 
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FIG. 13: The two functions of ( p^ ) are plotted for < i^/tt < 5. The 
pump parameter 9 has the value of 47r. The intersections marked with a 
little circle o determine the maxima in the photon number distribution, the 
other intersections determine the minima. 



in the photon number distribution. The conditions (|6T|) acquire a more uni- 
versal appearance when we introduce a new variable d = fy/n. Then we 
have 

sin^('i9) > (J) joy between minima and maxima, ^^^^ 
sin^(^9) < (J) joy between maxima and minima. 



where the pump parameter 9 = (pyr/A appears naturally. This is another 
reason for using 6 rather than (p. In passing we mention that the critical "i? 
values, for which sin^(^9) = holds, play a great role in the semiclassical 

theory of the micromaser. In the context of the quantum treatment discussed 
here, they determine the n values of the extrema in the photon number 



distribution. This is illustrated in Fig. [T^. The corresponding histogram of 
p{n) should show four maxima, and in Fig. 14 it does indeed. Of course, not 
all maxima have to be as pronounced as the two in the 6 = 2.097r case of 
Fig. |ll|(d). Again we note that, as in (|60D, the thermal photon number u is 
irrelevant in (|62D . The number of the extrema in the histogram of the photon 
number distribution as well as their locations are solely determined by the 
value of 6. Naturally, further details may depend strongly on the value of u. 

In the discussion of the plots in Fig. |] we had remarked that the pattern 
roughly repeats itself with a 6 period of 27r. That may be true for the 
parameter range of Fig. 0, but for other values of the effective rate r/A, other 
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FIG. 14: Histogram for a steady state with 6 = Att (actually: r/A = 36, u 



0.15, (p = 27r/3). The dashed line shows lOOp(n). Consistent with Fig. |T3 
we note four maxima: p(l) = 0.00020, p(7) = 0.11316, p(17) = 0.19397, and 
p(29) = 0.00029. 

thermal photon numbers u, or other 6 ranges matters tend to be different. As 



an example we take a look at Fig. p!5| where the rate is much smaller, namely 
r/A = 6, and the temperature much higher, z/ = 2. Rather than the regular 
pattern of Fig. ^ we observe an irregular behavior. The initial oscillation is 
quickly damped away, but builds up again at larger values of the Rabi angle 
if, that is: at longer interaction times. This phenomenon is reminiscent 
of and related to the so-called Jaynes-Cummings revivals. Strictly speaking, 
the latter occur when a two-level atom is in permanent (resonant) interaction 
with photons of one field mode. The Rabi angle is then proportional to the 
elapsed interaction time. Although the 0AM plots in Fig. [T^ refer to many 
atoms with a fixed common interaction time, it is nevertheless clear that the 
revivals in this figure are similar — both in their appearance and in their 
physical origin — to those Jaynes-Cummings revivals. 
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FIG. 15: Mean photon number and Fano factor in the 0AM steady state 
for an effective rate oir/A—Q) and a thermal photon number oiv — 2, as a 
function of the Rabi angle (/? in the range from (/? = to </? = Stt. 
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One-atom maser: Steady state. Energy balance 

One cannot measure the photon state itself in micromaser experiments — 
by counting photons, say — because that would violate the requirement of 
a very large quality factor of the resonator. The experimenter can only look 
at the atoms emerging from the resonator after having interacted with the 
photon field. The experimental data consists of the statistics of the final 
atomic states. For instance, one could count which fraction ends up in the 
excited — state and which fraction emits a photon and comes out in the 
deexcited ^: state. It is an essential task of 0AM theory to compute these 
fractions, and many other statistical properties of the emerging atoms, and 
to relate them, whenever possible, to the properties of the photon state inside 
the resonator. 

The calculation of the probability for a — ^ transition during the pas- 
sage of the atom through the resonator is based on the result (|25|). We recall 
that \—,n) — >■ |m, n + 1) happens with probability sin^i^ipy/n + 1). Therefore, 
if the atom encounters a state 

oo 

p[a^a) = \n)p{n) {n\ , (63) 

n=0 

then the total transition probability is given by 

oo 

E p{n) sin {(fy/n + 1) 

n=0 

tr{psin2[(^(aat)i/2]} (64) 
(sin2[(^(aat)i/2]) . 

This probability can be computed immediately once the steady state p has 
been determined in accordance with (^61). In a series of challenging experi- 
ments, actual measurements have been performed (see and the references 
therein, in particular and [1^), and very good agreement with the theo- 
retical predictions has been found. 

When we multiply the transition probability ([6^ ) by the beam rate r 
and the energy per photon hu, we obtain the rate at which the photon field 
gains energy: rhuo proh{— —>■ zii). In view of Eq. (|50|), the energy loss rate 
is likewise given by the product of the decay rate A, the energy per photon 
huj, and the difference between the actual mean photon number (a^a) and 
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the thermal photon number v: Ahuj{{a^a) — z^). In steady state these two 
rates must be equal. After canceling the common factor of huj we thus arrive 
at a fifth central result of micromaser theory: 



In the steady state of a one-atom maser, 
pumped by resonant — atoms, the energy 
balance implies the equality 

r{sin^[ip{aa^f'^]) = A{{a^a)-u) . 

V ' V ' 

oc energy gain rate oc energy loss rate 

The expectation value on the left-hand side 
can be measured experimentally. 



We emphasize that by determining the fraction of atoms that have emitted a 
photon one can experimentally measure the mean photon number inside the 
resonator. The link between the two quantities is provided by the theory of 
the 0AM and stated in (H). 

The way in which we derived this statement simply appealed to the phys- 
ical significance of the expressions on the two sides. A more formal derivation 
employs the two-term recurrence relation (^) in the form 

r sin^[if{aa}f/'^]p{a^a) = A[{u + l)p{a^a + 1) - u pia'' a)]{a'' a + 1) . (66) 

Upon taking the trace on both sides, we find 

r{sm^[if{aa^y/^]) = A{u + l){a^a) - Au{a^a + 1) , (67) 

which, after an elementary simplification, is recognized to be identical to the 
statement in (|65|). 

One-atom maser: Trapped states 

At zero temperature there are no thermal photons: z/ = 0. Then the contri- 
butions [c] and [e] are not present in Eqs. (|53| ) and (|^) and in Figs. and |[ 



(65) 
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FIG. 16: Histograms for the 0AM steady states obtained for r/A = 10, 
if = tt/2, and (a) u = 0, (h) u = 0.1. In case (a) the trapping limits the 
possible photon numbers to the values n = 0,1, 2, 3. No trapping occurs in 
case (b). 



As a consequence, the steady state ( ^q) of the one-atom maser is particularly 
simple, viz. 



p(ss)(ata) 1 = p(ata) = p(0) J] sin^(^v^) • 



[a'a 



(68) 



n=l 



In contrast to the general u ^ case of (|56D , here the individual factors of 
the product might vanish. Thus we note this property of (|68|) : 



(69) 



if sin'^i'^VN) = holds for = 1, 2, 3, . . . 
then p{a^a) = holds for a^a > A^ . 



One says that the photon number is 'trapped' below a^a = A^. 

Such trapped states are very characteristic of the 0AM dynamics. For 
purely energetic reasons, the mean photon number (a^a) could be as large 
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FIG. 17: Mean photon number in the 0AM steady states obtained for an 
effective pump rate of r/A = 10 at zero temperature (z/ = 0), as a function 
of the Rabi angle in the range < (/? < 27r. Please note the dips at = vr 
and if = 2'K that indicate trapped vacuum states. 

as the effective pump rate r/A [recall Eq. (]65|), here for i/ = 0]. The specific 
0AM dynamics, however, may limit the photon number to much smaller 
values. We illustrate this in Fig. |TB|, where the histogram for Lp = 7r/2, that 
is: = 4 in (^), shows nonzero probabilities only for the photon numbers 
n = 0,1,2,3, although the effective rate, r/A = 10, is more than twice as 
large. A second plot in this figure demonstrates how a small number of 
thermal photons (z/ = 0.1) removes the trapping. 

A very remarkable trapped state is the trapped vacuum state, which is 
realized if ip equals an integer multiple of vr. In Fig. |l^, the trapped vacuum 
ip = TT and cp = 2iT are visible as pronounced dips in the plot of the 
mean photon number. Note also the minimum at = vr/ that corresponds 
to a trapping below a^a = 2. 

When the vacuum is trapped, the atoms perform one or more complete 
Rabi cycles while traversing the resonator. In other words, each atom first 
emits a photon and then reabsorbs it — perhaps repeatedly — eventually 
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leaving the resonator in the excited state — in which it entered. This is, 
of course, the n = situation of Eq. (p7[) . Under these circumstances it 
does not matter how frequent the one-atom events are, nor if they occur in 
accordance with the Poisson statistics of (^). It is simply impossible to get 
photons into the resonator. 



Approximations in one-atom-maser theory 

A direct experimental demonstration of the trapped vacuum state has not 
been achieved as yet. To understand why, we recall that the theoretical 
result (^) is based upon various approximations. Let us reconsider them 
and their significance for the trapped vacuum state. In this context, the 
most important approximations are: 

The interaction time is assumed to be the same for all atoms. 
This assumption has two aspects: 
The atoms should have the same velocity v. 
Differing velocities lead to differing Rabi angle ip. Fortunately, however, 
the trapping condition requires ip = n, where sin^((^) has a minimum. 
Therefore, deviations from the ideal (p value of vr do not change anything 
to first order. A velocity control to about one percent is good enough, 
and that can be realized experimentally. 

Atomic decays to other levels do not occur. 
The Rydberg states of Fig. |^ do not meet this condition sufficiently well; 
a few percent of the atoms make ultraviolet transitions to low lying levels 
and thus stop interacting with the photons in the privileged mode before 
they actually leave the resonator. Therefore, one has to use Rydberg 
states that live even longer, the so-called circular states in which both 



the angular and the magnetic quantum number are maximal [0. The 
masing of such circular states has already been demonstrated experimen- 
tally [|12|. A one-atom-maser operation is likely to be realized in the near 
future. 

The transition frequency Q of Eq. (P) is assumed to be equal to the photon 
frequency u for all atoms. 

This has also two main aspects: 
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The resonator must he mechanically stable. 
Apparently one does not have to worry much about this. I am told that 
the detuning resulting from mechanical instabilities is less than 100 Hz.0 

Electric stray fields must he small. 
Here the concern is about Stark shifts of the atomic levels. Inasmuch 
as the stray fields are prominent only in the entry and exit ports of 
the resonator, and not so much inside where the relevant atom-photon 
interaction takes place, they are not very critical. In other experiments, 
however, in which the atoms are prepared in coherent superpositions of 
— and stray fields are a great nuisance because the differential Stark 
shift changes the phase relations between — and ^: and thus may destroy 
the desired superposition. 

The photon damping is neglected during the atom-photon interaction. 
This is harmless. Indeed, the Jaynes-Cummings model can be extended 



by the inclusion of photon damping. One then finds that the trapping 



condition of (^) is altered slightly, but trapping still exists. 

Collective events are not taken into account. 
This is a critical point. We shall deal with it in some detail. 

Two-atom events 

A typical number for AL/v in the Garching micromaser experiments is 10~^7r. 



For an effective beam rate oi r/ A = 10, the probability (|30D for a one-atom 
event is then exp(— 2rL/t') = exp[— 2(r/A) x [AL/v)] = 99.4%. In this 
situation, only 0.6% of the atoms participate in collective events, and most 
of those in two-atom events. 

Although these collective events are very rare, they may cause dramatic 
effects nevertheless. In particular, they destroy the trapped vacuum state 
completely. The physical reason is the following. After two atoms, say, have 
been in the resonator simultaneously, the photon field is no longer in the 
vacuum state, there are nonzero probabilities for one or two photons. Single 
atoms that follow will then emit more photons and so drive the maser field 



detuning should not be compared with the transition frequency of 21.5 GHz, but 
rather with the Rabi frequency of ~ 14 kHz. 
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FIG. 18: Histograms of (a) the trapped vacuum state; (b) the state im- 
mediately after the passage of an atom pair; and after subsequent one-atom 
operation for a period during which (c) two, (d) five, (e) one hundred, or (f ) 
two hundred one-atom events occur on average. 



further away from the vacuum. This is demonstrated in Fig. ^ Eventually, 
after a large number of one-atom events, the maser field will return to the 
vacuum state. The analysis given below shows that very many one-atom 
events are needed to undo the damage caused by a single two- atom event, 
so many indeed that another two-atom event is very likely to occur in the 
meantime. Therefore, we have to conclude that, despite their rare occurrence, 
the collective events prevent the observation of the trapped vacuum state in 
a micromaser with a poissonian pump beam. 

Consider the 0AM dynamics of Eq. (|53|) at zero temperature {u = 0) and 
for if = n, 



Q^pt[a'a) =ri sm 

— |ri sm 



n[a'a 

2 



Pt{a^a - 1) + A{a^a + l)pt{a^a + 1) 



(70) 
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where ri is the rate of one-atom events; its relation to the beam rate r is 
given by ri = r exp{—2rL/v). One verifies easily that the vacuum state 
Pvac = S{a^a, 0) is the steady state of ([7D|). 

For (f = n, the trapping condition is not only obeyed by = 1, but 
also by = 4, 9, 16, . . ., so that the threshold from three to four photons 
cannot be crossed by one-atom events. Now, after a two-atom event has 
happened to an initial vacuum state, the probabilities for one or two photons 
will be nonzero, but those for three or more photons still vanish. As long as 
only one-atom events occur thereafter, the probabilities for more than three 
photons remain zero. Thus a fitting ansatz for pt{a^a) is 

= ^5(ata, 0) + ap^(5(ata, 1) + ai^'>6{a^a, 2) + af^'>6{a^a, 3) , (71) 



wherein the time-dependent numerical coefficients a 



(0,1,2,3) 



are the respective 



probabilities to have no, one, two, or three photons in the resonator. The 
sum of these probabilities is 100% — this is the normalization of p to unit 
trace — and that requires = 1 — a^^^ — a^^-* — a^^\ 

(3) 

At t = 0, we have Og = 0, of course, and typical numbers for the initial 



values of al^'^ and al'^' are [|4 



„(2) 



sm 



ttVq)] = 0.3250, a^o 



(2) 



9 



sin(7rJ3/2) 



0.1575. (72) 



As implied by ([70|), the time dependence of these coefficients is determined 
by 



Adi 



( l + ^sin^l 
A 



v^) 







V 



■^sin2(7rv^) 
A 





2 + ^sin2(7ry3) 
A 

— ^ sin^(7r-\/3) 



I 



( .,(1) \ 



a 



(2) 



V 



a 



(3) 



0. 



(73) 



The time scale for the return to the vacuum is set by the eigenvalues of this 
3x3 matrix, more precisely: by the smallest one. For r/A = 10, ri/A = 9.937 
these eigenvalues are 



13.974, 6.7225, 0.063869 



ri/A 
155.6 



(74) 
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FIG. 19: Time dependence of the emission probability of atoms in one- atom 
events after the passage of an atom pairs, for < t < 400/ri. The abscissa 
is hnear in (ritY^^. 



Please note that the smallest one is very small indeed. It takes hundreds of 
one-atom events to get close to the vacuum state; see Fig. ^ again, where 
the photon states are reported for rit = 0, 2, 5, 100, 200. During such a long 
time, another collective event is very likely to occur, so that the vacuum state 
will not be reached at all. 

To be more specific, we note that the experimental signature for the 
trapped vacuum state is the observation of all atoms emerging in the excited 

state. Therefore, an appropriate numerical criterion should involve the 
probability that the atom emits a photon if it traverses the resonator at 
time t. This probability is given in Eq. (|63). For states of the form (|7l|) it 
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reads 

prob,(:!: ^ = ^ sm2(7rV2) + af^ sin^nVS) . (75) 

This probability is plotted in Fig. |1^. If we now require that 99% of the 
atoms emerge in the upper state —, then the emission probability must get 
as small as 1% before we would say that the vacuum state has been reached 
again. The time t = 370/ri has to go by for the parameters we have been 
using consistently, and the probability that at least one more collective event 
takes place in the meantime is 68.8% [Q. This number is discouragingly 
large — it is clear that one needs to operate the micromaser with a pump 
beam that has considerably fewer collective events, if the trapped vacuum 
state is to be observed. 

The small eigenvalue of the 3x3 matrix in ([73| ) is one theoretical indica- 
tion that very long time scales are present in the dynamics of the micromaser. 



In this context, it is worth mentioning that very recent experiments [15] have 
demonstrated the presence of extremely large time scales in rather different 
dynamical regimes of the micromaser as well. 

Correlations among emerging atoms 

We have remarked above that the photon state cannot be observed directly in 
micromaser experiments. The only reproducible data concerns the emerging 
atoms, more specifically: their statistical properties. Experimental tests of 
the micromaser theory must therefore compare the measured atom statistics 
to the theoretically predicted ones. We have seen one elementary example 
already, in the context of Eq. (p5|). On the one hand, one can calculate 
the transition probability (^); on the other hand, it can be measured by 
counting the atoms that end up in the — or the ^: state. 

More subtle are statistical predictions about correlations among emerging 



atoms. The general theory for their calculation is rather involved [|T6[ and well 
beyond the scope of this tutorial review. We shall therefore be content with 
presenting two rather simple examples. They serve the purpose of conveying 
the spirit in which such calculations are done without requiring too much 
technical detail. 

Consider a 0AM at zero temperature {u = 0) operated with a Rabi angle 
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statistical predictions, if nothing else is known. For instance, the (a priori) 
probability for an atom to end up in the ^: state is 



as stated in Eq. ([6^). 

Now suppose that we have just registered, at t = 0, an atom in the ^: 
state. What is then the probability that the next atom emerges also in the 
:^ state? We are here asking for a conditional probability, rather than for 
an a priori probability. The steady state p^^^^ is not appropriate any more. 
In its stead we must employ that time dependent state pt which correctly 
accounts for the imposed conditions. Then 



is the transition probability if the next atom comes after the elapse of time t. 
And the probability that the next atom arrives between t and t + dt is 
dtrexp{—rt), where we remind ourselves of the Poisson formula (pS]) . Upon 
putting things together, we find that the asked-for probability that the next 
atom also emits a photon is given by 



We still have to say how the correct pt is calculated. This requires to find 
(i) the appropriate initial state po and (ii) the appropriate master equation 
for Pt. 

At t = 0, we know that an atom has just been registered in the ^: state. 
Prior to its passage through the resonator we had no additional information, 
so that then p^^^^ applies. Immediately after the observation of the initial ^: 
atom, however, the photon state is reduced to account for this information. 




(77) 




(78) 




(79) 



state reduction: p' 



,(SS) 



6{a^a, 1) = Po • 



(80) 
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In other words: after the registration of the atom emerging in the ^: state 
we know with certainty that there is one photon in the resonator. 

Until the next atom arrives, there are no other atoms traversing the res- 
onator. Consequently, pt is determined by the master equation ( ^6]) of the 
unpumped resonator (here for z/ = 0). So we find first 



= (l - e-^') 5{a)a, 0) + e-^'5{a)a, 1) , (81) 

then 

prob,(- ^ ^:) = (l - e~^') sin2(7r/V2) , (82) 

and finally 

prob(next in zn) = ~ ^ sin2(7r/V2) . (83) 

This answers the question asked above. Compared with the a priori transition 
probability in (|77|), this conditional transition probability is smaller. 

The time dependent transition probability (|8^) vanishes for t = and 
approaches sin^(7r/ a/2) for t — oo. Both properties are easily understood on 
physical grounds. At t = 0, that is: immediately after the first atom was 
found in the state, the trapping condition forbids the emission of another 
photon into the resonator. And if no atom has traversed the resonator for 
a time long compared to the photon life time 1/A — this is the physical 
significance of the mathematical t oo limit — the photon field is in the 
vacuum state, and then the transition probability is given by ( pSD with n = 
and if = T[ / 1/2. 

As another example, we consider again that at t = an atom is found in 
the ^: state, but let us now ask this question: If another atom comes between 
t and t + dt, with which probability will it be in the ^: state as well? Note 
the difference. Above we were interested in the next atom, now in another 
atom. Clearly, the initial state is obtained by the same reduction, viz. that 
of (|80D , but now pt is determined by the 0AM master equation (|70D with 
= Tx I \f2. Rather than (p2D, we get here 



Asin^(7r/V2) . 

prob,(=!: -^^) = -— .2. I fK. (1 - 

yl + r sm (7r/\/2) ^ 



- L4 + rsin^('7r/V2) U 



a priori probability correlation function for m atoms 

(84) 
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We recognize that the first factor is the a priori probabihty of ([77|), and the 
second factor is the so-called correlation function. It is a numerical measure 
for the strength of the correlations among atoms emerging in the zi: state. 
If the temporal separation is too large — physically: t ^ l/A,l/r — there 
must not be any correlations. Indeed, the correlation function approaches 
unity for t — > oo, and the conditional probability (Q) equals the a priori 
probability in this limit, as it should. 

These two examples must suffice as illustrations for the typical reasonings 
that are necessary for the computation of statistical properties of the emerg- 
ing atoms. Let us just mention that the experimental data refers actually 
to the statistics of the detector clicks, rather than to the atoms themselves, 
and therefore the efficiencies of the detectors have to be taken into account 
as well. Naturally, that introduces further complications which are, however, 
under good theoretical control. 
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